GENERAL RELATIVISTIC EFFECTS ON 
THE INDUCED ELECTRIC FIELD 
EXTERIOR TO PULSARS 

S. SENGUPTA * 

Indian Institute Of Astrophysics, Bangalore 560 034, India 
To appear in The Astrophysical Journal, August 10, 1995 



*e-mail : sujan@iiap.ernet. in 



1 



Abstract 

The importance of general relativity to the induced electric field exterior to pul- 
sars has been investigated by assuming aligned vacuum and non-vacuum mag- 
netosphere models. For this purpose the stationary and axisymmetric vector 
potential in Schwarzschild geometry has been considered and the correspond- 
ing expressions for the induced electric field due to the rotation of the magnetic 
dipole have been derived for both vacuum and non-vacuum conditions. Due to 
the change in the magnetic dipole field in curved spacetime the induced electric 
field also changes its magnitude and direction and increases significantly near the 
surface of the star. As a consequence the surface charge density, the acceleration 
of charged particles in vacuum magnetospheres and the space charge density in 
non- vacuum magnetosphere greatly increase near the surface of the star. The re- 
sults provide the most general feature of the important role played by gravitation 
and could have several potentially important implications for the production of 
high-energy radiation from pulsars. 

Subject headings: electromagnetism : theory - pulsars : general - relativity - 
stars : neutron 
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1. INTRODUCTION 

The realization that pulsars are most likely to be neutron stars surrounded by a dense 
magnetosphere (Gold 1968; Pacini 1968) and the pioneering work by Goldreich and Julian 
(1969) have provided the basic physical model on which virtually all particular models have 
been based. At present two general types of models are being considered in-order to un- 
derstand the phenomenon of high energy radiation from pulsars viz., the polar cap model 
(Arons 1983) and the outer gap model (Cheng, Ho & Ruderman 1986). Goldreich and Julian 
first realized that inspite of the intense surface gravity, the neutron star must possess a dense 
magnetosphere. Unfortunately the actual role played by gravitation in this aspect has not 
been investigated so far in detail. 

General relativistic effects near pulsars have been explored only in the context of emission 
models for X-ray pulsars and gamma-ray burst. The effect of light bending in X-ray burst 
and X-ray pulsars was studied by Meszaros and Riffert (1988), Riffert, Meszaros & Bagoly 
(1989). Recently Gonthier and Harding (1994) have examined the importance of general 
relativistic corrections to the production of gamma rays near the surface of neutron stars. 

It is well-known that the induced electric field due to the rapid rotation of the magnetic 
dipole of a neutron star plays the main role in determining the characteristic of the pulsar 
magnetosphere since this quantity determines the number density of charged particles torn off 
from the surface of the star, the acceleration of charged particles and hence the production 
of electromagnetic radiation, etc. Since the magnetic field in curved spacetime differs in 
both magnitude and direction from that in flat spacetime, the induced electric field is also 
different in curved spacetime. Therefore it is very important to include general relativity in 
any discussion of pulsar magnetosphere. 

In this paper the importance of general relativistic corrections to the induced electric field 
external to a neutron star has been investigated quantitatively by considering the simplest 
aligned vacuum and non-vacuum magnetosphere models. 

The paper is organized in the following way: In section 2 the Schwarzschild metric and the 
corresponding orthonormal tetrad components alongwith the transformation rules are given. 
In section 3 we discuss the dipole magnetic field in Schwarzschild spacetime which has been 
used to determine the induced electric field in section 4. In section 4 a detail quantitative 
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discussion on the effects of spacetime curvature to the various important quantities governed 
by the electric field has been presented. For simplicity we have assumed the vector potential 
to be stationary and axisymmetric. Finally we conclude the discussions in section 5. 

2. THE SPACETIME METRIC 

Since we are interested on the electromagnetic field exterior to the neutron star, we con- 
sider the Schwarzschild metric as it provides a simple description of the spacetime geometry 
exterior to neutron stars. 

The metric is given by : 

ds 2 = (1 - 2m/r)c 2 dt 2 - (1 - 2m/r)- 1 dr 2 - r 2 d6 2 - r 2 sin 2 6d(f) 2 (1) 

where m = MG/c 2 , M being the total mass of the star. 

The non-zero components of the orthonormal tetrad of the local Lorentz frame for 
Schwarzschild geometry are given by 

X\ t) = (1 - 2m/r)- 1 / 2 , X\ r) = (1 - 2m/ r)^ \\ e) = \, = 

If -F( Q /3) are the components of the electromagnetic field tensor in local Lorentz frame, 
then the components of the electromagnetic field tensor Fy defined in the curved spacetime 
through 

F(af3) = \\ a) \\ 0) F i:j (2) 

3. THE DIPOLE MAGNETIC FIELD IN CURVED SPACETIME 

Since the magnetic field of a neutron star is generally approximated to be that of a pure 
magnetic dipole, in flat spacetime we can write the components of the magnetic field (in 
polar coordinates) as 




where \x is the magnetic dipole moment. 

The magnetic dipole has to be modified in curved spacetime. In Schwarzschild spacetime 
the electromagnetic potential A { is given by (Petterson 1974; Wasserman & Shapiro 1983) 

Ai = (0,0,^,0) 
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with 

As = - 3/i Sm „ ® [r 2 ln(l - 2m/ r) + 2mr + 2m 2 ] (4) 

Using the definition = {Aj ti — A it j) we get the components of the magnetic field in 
Schwarzschild geometry as 

F r4> = - 3/ " Sm 9 2g [- ln(l - 2m/r) + (1 - 2m/r)- 1 + 1] (5) 
Am z m 

sin#cos#[ln(l — 2m/ r) H (1 + m/r)] (6) 



4m 3 r 
Therefore using equation (2) we get 

BJSch) = - 3/ " C °^ [ln(l - 2m/r) + — (1 + m/r)} (7) 
4m d r 

Sef^c/i) = 3/iS1 9 ng [- ln(l - 2m/r) + (1 - 2m/r)- 1 + 1](1 - 2m/rf/ 2 (8) 

It should be mentioned that the above expressions are specialized to the lowest order 
expansion in terms of spherical harmonics. Nevertheless they are valid in the context of 
neutron stars since the magnetic field is not so strong relative to the gravitational field such 
that the spacetime geometry could be affected significantly by the magnetic field. 

4. THE INDUCED ELECTRIC FIELD IN CURVED SPACETIME 

In order to investigate the effect of curved spacetime on the electric field exterior to a 
rotating neutron star, we consider two simplest models, viz., (1) the aligned vacuum rotator 
and (2) the aligned non- vacuum rotator. Since the effect of curvature would be the same, the 
results will provide the most general feature for any magnatospheric model which assumes 
stationary (J^ = 0) and axisymmetric (A = 0) behavior of the vector potential. In the 
present study we have ignored any effects due to magnetosphere plasma. Also we have 
neglected the dragging of inertial frame due to rotation which is insignificant at the exterior 
of even millisecond pulsars. 

4.1. The Aligned Vacuum Rotator 

Since the derivations are nothing but a generalization of Goldreich - Julian solution, we 
briefly discuss it before going to derive the electromagnetic field in curved spacetime although 
it can be found in any text book on Pulsars e.g., Shapiro and Teukolsky (1983). 
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First of all it is considered that a rotating neutron star has an aligned dipole external 
magnetic field given by equation (3). The stellar material is assumed to be an excellent 
conductor so that just inside the star an electric field will be present which satisfies 

E int + < ^2il x B int = (9) 
c 

where uj is the angular velocity vector of the star. 

Here and afterwards the superscript 'int' is to be understood to represent the quantities 
just inside the star and not far bellow the surface where the above condition may not be 
valid. 

Assuming there are no surface currents, both the normal and the tangential components 
of B are continuous across the stellar surface. Thus just inside the surface the magnetic field 
can be written by taking r = R in equation (3) where R is the radius of the star. Then 
equation (9) gives for E just inside the surface 

E in *= sin 2 6,- sine cos 6,0) (10) 

cR 2 2 

The tangential component of E is continuous across the surface, so just outside the star 
equation (10) implies 

where P<i is the Legendre polynomial of second degree. 

In the present case since the exterior is assumed to be vacuum, so 

E ext = 

where 

V 2 = (12) 

In order to satisfy the boundary condition given by equation (11) at r — R, the solution 
of equation (12) must be 

2uuR 2 „ . 



6 



and hence 

W xt = -^^[(3cos 2 £-l),2sintfcos£,0] (13) 

Now we will derive the electric components exterior to a rotating neutron star for 
Schwarzschild background geometry with the basic assumptions taken in the above dis- 
cussions unaltered. 

Equation (9) can be written in its covarient form as 



F s p u d = (14) 



where u s are the components of the four velocity vector. 
From equation (14) we obtain 



u v 

F H = F^ r — (15) 



Fet = (16) 



(Note that u = (u*, 0, 0, u^) for the present case.) 

Using equations (5) and (6) and noting that in Schwarzschild geometry 



E r — F( rt ) — F rt 



Eg = F m = F et /[r(l - 2m/r) 1 ' 2 } 

and taking 

to 

u l c 

we obtain the induced electric field components just inside the surface 
(r ~ R) as 

Ei = 3/^sin 2 g _ + L _ + 

= 3^ cos g sing _ + 2m + m _ 2m 1/2 

4cm ,i " r r r 
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If f(r) and g(r) are two arbitrary functions which solely determine the spacetime curva- 
ture effect then in curved spacetime we can write 



Ei nt (curved) = ^ sin 2 9 fir) (19) 
cr 



Ei nt (curved) = sin9 cos 9 g(r) (20) 



Comparing equations (19) and (20) with equations (17) and (18) we obtain 

f(r) = |4 [(1 - 2m/r)~ 1 + - ln(l - 2m/ r) + 1] (21) 
Am 2 m 

^ r = -5-3NI - 2m r) + — 1 + - 1 - 1/2 22 

8m d r r r 

In the above derivations we have assumed that the metric can well describe the spacetime 
curvature just bellow the surface (r ~ R) of the star, i.e., the curvature effect just bellow 
the surface is the same as that at the surface to the star. 

In flat spacetime the components of the exterior electric field is given by equation (13). 

Therefore in Schwarzschild geometry the components of the exterior electric field can be 
written as 



E e r xt (Sch) = -^-(3cos 2 #- l)f(r) 



A ' !j;R ; (3cos 2 9- 1)[(1 - 2m/r)~ 1 + — ln(l - 2m/r) + 1] (23) 



Acm 2 r 2 m 



E e e xt (Sch) = _ 2 ^ 2 sin cos g (r) 



cr A 

?2 



3fiuR cos 6 1 sin ^ . . 2m m \i/i ^ m 



. , [ln(l - 2m/r) + — (1 + -)](1 V' 2 (24) 

4cm 3 r r r r 

Clearly E ext (curved) — > E ext (flat) as r — > 00 as desired. 

Taking as examples, the Crab pulsar with mass M equal to 1.4 solar mass, radius R equal 
to 10 6 cm, period P equal to 33 ms and the magnetic dipole moment \i equal to 2 x 10 30 Gcm 3 
we show in Figure 1 the r component of the quadrapole electric field at an angle 9 = 0° 
to the axis of rotation as a function of distance from the stellar surface for flat and curved 
spacetimes. It is found that a significant increase in the electric field intensity near the 
surface is caused due to the inclusion of general relativistic effect. 
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Now we have the functional form for the dipole magnetic field and the quadrapole electric 
field in flat spacetime as well as in Schwarzschild spacetime. We would like to see how the 
field lines differ in two different geometries. The field lines represent a constant magnetic 
and electric flux b{9) and e{9) respectively and can be found by rotating a given field line 
about the z-axis in a spherical co-ordinate system with the magnetic dipole aligned with the 
z-axis. The fluxes from 9 = to 9 are given by 

6(0) = f B r .dsi= f B r r 2 dtt (25) 
Jo Jo 

e(9) = [ E r ..da= [ E r r 2 dQ (26) 
Jo Jo 

The magnetic field lines in flat spacetime and in curved spacetime have been presented 
by Gonthier and Harding (1994). 

For an induced quadrapole electric field in flat spacetime the field lines are 

lniiujR 2 . , . . 

e{9) = — sin 9 cos 9 = constant (27) 

cr 2 

In curved spacetime the field lines are given by 

= _ SirpuR 2 s . n2 Q _ 2 ! + r_ _ + 

2m 2 c m 
= constant (28) 

Taking the Crab pulsar as an example with a period of 33 ms we show in Figure 2a the 
quadrapole electric field lines in the z-x plane out to 40 stellar radii, beyond which the effect 
of curvature becomes insignificant. Figure 2b shows a closer view of the same field lines. 

Since the normal component of the electric field is discontinuous for the case of aligned 
vacuum rotator, the corresponding surface charge density can be written in flat spacetime 
as 

p = ±.(ET - E^) = cos 2 9 (29) 
In curved spacetime the same quantity can be written as 

3fIU cos 2 9[(l-2m/R)~ 1 + - ln(l - 2m/i?) + 1] (30) 



8ircm 2 m 
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Since the effect of curvature increases the electric field intensity near the surface, therefore 
in curved spacetime the surface charge density is much higher than that in flat spacetime. 
For a pulsar with mass 1.4M Q and radius 10 km the surface charge density at the pole in 
curved spacetime is twice of that in flat spacetime. 

The Lorentz invariant scalar product E ext .B does not vanish. This quantity gives a 
measure of the force which a co-rotating charged particle feels in the direction of the magnetic 
field. 

In flat spacetime 

E-.B = _V^ Wg (31) 

leading to an acceleration of a charged particle in the direction of the magnetic field: 

eE ext .B A - "" d2 



a = 



m p |B| rripC r 4 



^_f^_ cog 3 (3 cog 2 6 + 1} -l/2 ( 32) 



where m p is the mass of the charged particle and e is its charge. The corresponding quantities 
in curved spacetime can be obtained by using equations (7), (8) and (24) and are given by 



9fi 2 uR 2 



E B = ^V F(r)G(r)c ° S 6 (33) 
and 



3e uuR 2 , o ^cos 2 6 ~, . sin 2 9 , 2m. ~, .,■,/-, . . 
^f(r)G(r) cos 3 8[ ^G 2 (r) + — =-(1 )F 2 (r)]- 1 / 2 (34) 



2m p c m 3 r 2 m 2 



where 



F(r) = [(1 - 2m/r)- 1 + — ln(l - 2m/r) + 1] (35) 

G(r) = ln(l - 2m/r) + — (1 + -) (36) 

r r 

In Figure 3 we present the variation of acceleration of charged particles with distance 
by considering the Crab pulsar with a rotational period 33ms. The acceleration has been 
measured at the angle 6 = 0° to the axis of rotation. Since gravitation greatly increases the 
electric field intensity near the surface of the star, the acceleration of charged particles is 
much higher in curved spacetime than that in flat spacetime near the surface. This result 
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indicates that the radiation spectrum of pulsars may significantly be influenced by the intense 
gravitational field of neutron stars and requires further attention. 

Since in realistic situation the magnetosphere must be filled up with charge plasma, there- 
fore the various quantities discussed in this section may not represent the actual situation. 
However the effect of curvature remains the same whether the region under consideration is 
vacuum or not, since the gravitational effect due to the magnetospheric plasma is negligible. 

4.2. The Aligned Non-vacuum Rotator 

It is now well accepted that the magnetosphere of neutron stars must be filled up with 
plasma such that E.B = everywhere within the light cylinder. In the Goldreich- Julian 
model the source of plasma is assumed to be free field emission from the surface. 

In this section we will briefly discuss the effect of spacetime curvature on the induced 
electric field by considering the aligned non-vacuum rotator model. 

In order to investigate the effect of curvature to the induced electric field we consider 
the Go dreich- Julian solution which simply corresponds to extending the interior solution to 
infinity with space charge density continuing smoothly through the surface. We neglect, for 
simplicity any effects due to the magnetospheric plasma. However inclusion of the same 
would not change the general feature of the curvature effect. 

Therefore, for this case, equations (17) and (18) express the components of the electric 
field outside the star in Schwarzschild spacetime with r being extended from just inside the 
surface to infinity (the superscript 'int' should now be ignored). 

The corresponding expression in flat spacetime can be written as 



Asymptotically as r — > oo, E(curved) — > E(flat) as it must be the boundary condition. 
From equation (10) we obtain the electric field lines in flat spacetime as 



E(flat) = ^(sin 2 0, -2 sin 9 cos 9, 0) 



(37) 



e{9) 



[cos' 



3 9 — 3 cos 9 + 2} = constant 



(38) 



3c 



From equation (17) we obtain the same in curved spacetime as 



e(9) = ^I-[ C o S 3 9-3cos9 + 2}[{l-2m/ry 1 + — ln(l - 2m/r) + 1] 




= constant 



(39) 
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Clearly the electric field lines for flat spacetime is independent of r whereas those in 
Schwarzschild background geometry depend on r only near the surface. At comparatively 
large distance the field lines for curved spacetime asymptotically become independent of 
distance and hence it is not possible to plot them at a large distance. Also, since the electric 
field is continuous from the interior to the exterior of the star for the non-vacuum case, the 
expressions are independent of the radius of the star. 

Taking, as earlier, the Crab pulsar with a period of 33ms we show in Figure 4 the induced 
electric field lines in curved spacetime in the x-z plane out to 4 stellar radii. Beyond this 
distance the field lines almost become independent of r and it is not possible to present them 
graphically. 

Figure 5 shows the complete picture of the non- vacuum magnetosphere in curved space- 
time with both the induced electric field lines and the dipole magnetic field lines obtained 
by using the expression for the magnetic field lines: 

, 3irur 2 sin 2 9 ri n , . 2m m,, . lnS 

b(9) = In 1 - 2m r) + 1 + — = constant 40 

4m 3 r r 

In both the plots we have taken, for convenience, the electric and the magnetic field lines 
from the interior of the star and extended it to the exterior. However, far bellow the surface 
the spacetime metric considered here, cannot describe the geometry of the region, hence in 
the plots they are not shown at a distance far bellow the surface of the star. 

From equations (27) and (28) we notice that the ratio between the electric fluxes in 
curved and flat spacetime for the vacuum model is the same to that for the non-vacuum case 
and is given by 

e(9)[curved} vac _ e(9)[curved} non - vac _ 3r 2 F(r) 
e(9)[flat] vac e(9)[flat] non _ vac 4m 2 

where F(r) is given by equation (35). 

Hence the effect of spacetime curvature to the electric field lines remains the same in this 

case as well. 

In flat spacetime the space charge density is given by 

a = —divE = ^- (3 cos 2 - 1) (42) 

47r 2ncr 3 

which is always quadrapolar. 
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The corresponding expression in curved spacetime can be written as: 



a 



16ncm 2 
2sin 2 1 



r 




2 e 




(1 - 2m/r)- 2 + 



(43) 



Figure 6 presents the space charge density as a function of distance from the surface of 
the Crab pulsar at an angle 9 = 0° to the axis of rotation. In all the figures the various 
quantities have been measured in CGS units and their absolute values have been taken in 
logarithmic scale unless mentioned otherwise. 

Since gravitation increases the electric field intensity at the surface, more particles get 
torn off the surface. The electric field required to pull ions from the molecular chains with 
the ions distributed in a one dimensional lattice along the chain and with an outer sheath 
of electrons, is of the order 1Q 12 V cm" 1 (Ruderman & Sutherland 1975). Therefore the effect 
of gravitation in some pulsars may cause the surface electric field to reach this cut off value 
and hence may enable to pull ions. However, ions may be torn off from the surface due to 
other effects (Mitchel 1982). In that case, the inclusion of gravitation would increase the 
number density of ions, in the same way it increases the number density of electrons. 



We have presented a quantitative discussion on the effects of spacetime curvature to 
the induced electric field of a pulsar. Considering stationary and axisymmetric electromag- 
netic vector potential in Schwarzschild background geometry the components of the induced 
electric field that arises due to the rotation of the magnetic dipole, have been derived and 
compared with that in flat spacetime. It is found that spacetime curvature increases the 
electric field intensity significantly near the surface of the star. As a consequence, in vacuum 
magnetosphere the acceleration of charged particles along the direction of the magnetic field 
increases significantly near the surface. Also the increase in the electric field intensity near 
the surface results in increasing the surface charge density significantly. The quadrapole 
electric field lines have been presented for Crab pulsar, as an example, in both curved and 
flat spacetime and it is found that the difference in the field lines for the curved and the flat 
spcetimes is well distinguishable upto 40 stellar radii. 

In a plasma filled magnetosphere the electric field lines are independent of distance in 



5. CONCLUSIONS 
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flat spacetime but in curved spacetime they depend on r near the surface. However, the field 
lines in curved spacetime asymptotically become independent of r at large distance. The 
induced electric field lines in curved spacetime have been presented by taking Crab pulsar, 
as an example and the combined dipole magnetic field lines and the induced electric field 
lines near the surface of the star for curved spacetime have been presented as well. 

Due to the increase in the electric field intensity near the surface the space charge density 
of a plasma filled magnetosphere increases significantly from that in flat spacetime. All 
the quantities measured in curved spacetime, however asymptotically reach their respective 
values in flat spacetime as a result of the usual boundary condition. Since the ratio of 
the components of the electric field in curved spacetime to that in flat spacetime is the 
same, therefore the effects of gravitation on various electromagnetic mechanisms remain the 
same for all types of magnetospheric models which assume stationary and axisymmetric 
vector potentials. Hence the present results are very important in the context of radiation 
mechanisms in pulsars. 

Though we have considered the most idealized models of the pulsar magnetosphere the 
effects of spacetime curvature to the induced electric field would remain the same provided 
there is no dipole or gravitational radiation. The effect of inertial frame dragging is negligible 
even for millisecond pulsars but the depurture from spherical symmetry due to rapid rotation 
of the star could alter the scenario significantly. Nevertheless, the present results would 
provide reasonably well insight on the effect of spacetime curvature to the induced electric 
field outside a rotating neutron star with its magnetic axis aligned to the rotational axis. 
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FIGURE CAPTIONS 



FIG 1. - The r component of the quadrapole electric field at an angle 6 = 0° to the 
axis of rotation for Crab pulsar as a function of distance from the stellar surface. Solid curve 
represents the field intensity for curved spacetime while dashed curve for flat spacetime. 

FIG 2a. — Quadrapole electric field lines (induced) for the Crab pulsar with fi = 
2 x 10 30 Gcm 3 . Solid curves represent field lines for curved spacetime while dashed curves 
for flat spacetime. The corresponding value assigned to the field line is the same for both 
curved and flat spacetimes. 

FIG 2b. — A closer view of the quadrapole electric field lines shown in Figure 2a. 

FIG 3. — The acceleration of an electron in the direction of the magnetic field at an 
angle 9 = 0° to the axis of rotation for Crab pulsar as a function of distance from the stellar 
surface. Solid curve - curved spacetime; dashed curve - flat spacetime. 

FIG 4. - The induced electric field lines in curved spacetime for the Crab pulsar with 
a period 33ms and fi = 2 x 10 30 Gcm 3 by considering the magnetosphere to be plasma filled 
(see text). 

FIG 5. — The combined dipole magnetic and induced electric field lines in curved 
spacetime for the Crab pulsar with /i — 2 x 10 30 Gcm 3 by considering plasma filled magne- 
tosphere. Solid curves represent the induced electric field lines while dashed curves for the 
dipole magnetic filed lines. 

FIG 6. — The space charge density at an angle 6 = 0° to the axis of rotation for the Crab 
pulsar as a function of distance from the stellar surface. Solid curve - curved spacetime; 
dashed curve - flat spacetime. 
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